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INTRODUCTION 

A  body  falling  through  the  atmosphere  may  be  used  to  measure 
the  density  of  the  atmosphere,  if  the  aerodynamic  characteristics 
of  the  body  are  well  determined  beforehand  and  the  velocity  history 
of  the  fall  is  measured  and  recorded.  A  spherical  shape  has  been 
utilized  to  probe  the  upper  atmosphere  of  the  earth,  since  it  has 
the  inherent  advantage  of  possessing  rotational  symmetry,  thus 
eliminating  the  problems  of  angle  of  attack. 

Newton* s  second  law,  written  for  a  falling  body  in  a  simplified 
form  (coriolis  acceleration,  buoyancy,  etc.,  neglected),  is 

2FV  =  mav 


where  v  denotes  the  velocity  direction.  Since  drag  is  the  only 
force  in  the  velocity  direction,  the  above  relation  becomes 

-Drag  =  mav 


Applying  the  definition  of  drag  coefficient 


Drag 


1 

2 


p  V2  A 

1 


one  gets 

p  V2  A  CD  =  mav 

which  can  be  solved  for  density 


\  V2  A  CD 

It  is  seen  that  the  density  of  the  surrounding  atmosphere  can  be  1 
obtained  if  values  for  the  terms  on  the  i*ight-hand  side  can  be 
measured.  The  mass,  m,  and  the  area.  A,  on  which  the  drag  coeffi¬ 
cient  is  defined  are  constant  and  known  before  the  flight.  There 
have  been  two  methods  used  to  determine  the  instantaneous  velocity, 

I 
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V,  and  the  acceleration  (or  rather  deceleration  in  this  case),  ay* 
in  the  direction  of  the  velocity.  The  acceleration  has  been  measured 
directly  using  built-in  accelerometers  and  then  the  velocity  obtained 
through  integration.  The  units  carrying,  these  accelerometers  have 
been  rather  heavy  and  have  successfully  been  used  at  lower  altitudes. 
To  obtain  measurements  at  high  altitudes,  much  lighter,  larger 
diameter  bodies  need  to  be  used  so  that  the  fall  rate,  as  well  as 
the  mass,  will  be  low  enough  to  detect  density  gradients  in  the 
tenuous  environment.  One  such  object  commonly  used  has  been  the 
balloon  with  a  metallic ; ref lector  inside  so  that  it  can  be  tracked 
by  radar.  With  precision,  ground-based  radar,  the  position  of  the 
falling  balloon  is  determined  at  each  instant  of  flight.  This  space- 
time  plot  is  then  differentiated  to  obtain  the  velocity  history  and, 
again* to  obtain  the  acceleration. 

One  more  quantity  needs  to  be  determined  before  density  can  be 
calculated.  !Rie  drag  coefficient  of  the  falling  shape  must  be  known 
at  each  instant  during  the  fall.  It  is  the  subject  of  this  report 
to  investigate  the  validity  of  the  available  data  on  spherical  shapes 
for  use  in  the  falling  sphere  technique  of  measuring  the  atmospheric 
density.  Spheres  have  been  widely  used  for  such  measurements,  and  a 
considerable  amount  of  experimentation  has  been  performed  on  this 
configuration  to  determine  its  drag  characteristics  in  different 
flight  regimes.  These  data  have  been  collected  and  are  analyzed  in 
this  report.  • 

In  continuum  flow,  the  drag  coefficient  of  a  sphere  has  been 
described  as  a  function  of  the  Reynolds  number  and  the  Mach  number. 

In  rarefied, gas  flow  regimes,  the  Knudsen  number  is  also  a  signifi¬ 
cant  parameter.  The  Knudsen  number,  however,  can  be  expressed 
approximately  as  a  Mach  number  -  Reynolds  number  function.  The 
ranges  of  these  parameters  that  are  of  interest  in  the  falling  sphere 
technique  of  measuring  the  upper  atmosphere  are  as  follows: 


1  Subsonic 

]  i  Mach  Number  =  less  than  one  . 

Reynolds  Number  =  5  x  1q2  to  5  x  10^ 

i  -  ,  ' 

Supersonic 

Mach  Number  =  1  to  5. 

:  Reynolds  Number  -  5  x  10^  to  5  x  10 ^ 

! 

These  regimes  are  illustrated  in  Figure  1  by  a  shaded  area.  The 
data  in  each  regime,  the  subsonic  and  the  supersonic,  will  be 
analyzed  separately. 

SUBSONIC  REGIME 

The  subsonic  data  may,  again,  be  subdivided  into  two  main 
categories,  the  incompressible  and  the  compressible.  The  incompressible 
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data  are  obtained  by  performing  experiments  in  an  essentially  incom¬ 
pressible  fluid,  like  water,  or  they  may  be  obtained  in  a  gas,  under 
conditions  such  that  the  compressibility  effects  are  negligible. 

The  upper  limit  under  which  a  gaseous  flow  customarily  has  been 
considered  incompressible  is  when  the  Mach  number  is  equal  to  or 
less  than  0.3.  Above  that  Mach  number,  the  drag  coefficient  is 
also  a  function  of  Mach  number  in  addition  to  being  a  function  of 
Reynolds  number  alone. 

Incompressible  Data 

Figure  2  shows  the  Reynolds  number  coverage  of  pasc  experiments 
available  in  open  literature  on  subsonic  sphere  drag  ana  the  range 
of  applicability  for  the  falling  sphere  atmospheric  sensing  technique. 
A  work  that  describes  the  drag  coefficient  over  a  wide  range  of 
Reynolds  numbers  is  that  of  Wieselsberger,  references  (1)  and  (2). 

In  this  investigation,  the  drag  coefficient  of  a  sphere  was  measured 
over  a  range  of  Reynolds  numbers  from  8  x  102  to  9  x  1CP.  These 
experiments  were  performed  in  a  wind  tunnel  at  Mach  numbers  less 
than  0.1  and,  together  with  Allen’s  measurements  in  water,  reference 
(3)j  at  Reynolds  numbers  below  2  x  102,  define  what  will  be  called 
the  Wieselsberger  curve  or  the  incompressible  curve. 

Transition  from  laminar  to  turbulent  flow  in  the  boundary  layer 
takes  place  at  Reynolds  numbers  approximately  between  2  x  105  an(j 
5  x  1CP,  which  is  accompanied  by  a  significant  and  a  rather  abrupt 
decrease  in  the  drag  coefficient.  This  phenomenon  has  attracted 
considerable  attention.  In  references  (4)  through  (8),  different 
phenomena  affecting  transition  in  the  boundary  layer  are  investigated. 
In  these  references,  the  effects  of  tunnel  turbulence,  Mach  number, 
surface  roughness,  etc.,  are  studied.  Tripping  of  the  boundary 
layer  and  its  effects  are  also  described,  and  the  resulting  changes 
in  the  pressure  distribution  around  the  body  are  measured  in  the 
critical  Reynolds  number  range.  Since  the  critical  Reynolds  number 
is  outside  the  range  of  interest  in  the  falling  sphere  technique 
of  atmospheric  sensing,  a  detailed  analysis  of  these  experiments 
will  not  be  made. 

The  range  of  Reynolds  numbers  of. interest  at  subsonic  Mach 
numbers  is  between  5  x  102  and  5  x  10^.  Seven  investigations 
have  been  performed  in  this  region,  satisfying  the  incompressibility 
condition,  and  they  are  described  in  references  (l),  (2),  (3),  {9)» 
(10),  (11),  and  (12).  All  of  the  measured  data  points  from  these 
references  are  reproduced  on  a  Cd  versus  R^  plot  in  Figure  3.  The 
solid  line  is  reproduced  from  Wieselsberger' s  paper,  reference  (1), 
as  well  as  the  data  points  substantiating  this  line  from  references 
(2)  and  (3).  Wieselsberger' s  experiments  were  performed  in  a  wind 
tunnel  with  the  model  supported  on  wires.  The  force  on  the  sphere 
was  determined  by  the  amount  the  model  is  displaced  on  the  support 
system  against  counterweights. 

The  other  incompressible  data  in  air  were  obtained  by  dropping 
spheres  from  towers  (Shakespeare,  reference  (12))  and  in  mine  shafts 
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(Lunnon*  reference  (10)).  These  data  differ  considerably*  especially 
the  ones  obtained  in  mine  shafts  from  the  wind-tunnel  experiments. 
This  can  be  attributed  to  the  difficulties  that  are  encountered  when 
timing  the  fall  of  a  body  over  a  certain  height.  As  is  pointed  out 
by  the  author  in  reference  (12)*  the  accuracy  of  the  experiments 
heavily  depended  on  the  atmospheric  conditions,  The  most  consistent 
results  from  the  tower  drop  tests  were  obtained  on  dull*  cloudy* 
and  still  days.  Clear  weather  invariably  brought  more  scatter  into 
the  data.  This  was  attributed  to  drafts.  The  models  used  were 
hollow  celluloid  spheres*  which  were  weighted  with  lead  to  obtain 
the  desired  mass.  Sphericity  of  these  spheres*  which  were  such 
that  they  could  be  opened  to  insert  the  lead*  can  also  be  questioned. 
The  data  obtained  in  mine  shafts  (Lunnon*  reference  (10))  are  belo w 
the  wind-tunnel  measurements  by  as  much  as  20  percent.  Here  again* 
air  currents  in  the  shaft*  as  well  as  temperature  variation  with 
depth*  could  be  contributing  to  the  inaccuracy  of  the  measurements. 

There  are  three  references  describing  results  of  measurements 
performed  in  water.  In  1900*  Allen*  reference  (3),  measured  drag 
by  dropping  small  steel  spheres  in  water.  The  measurements  are 
shown  with  a  symbol*^*  in  Figure  3  and  are  seen  to  be  from  10 
percent  to  17  percent  below  the  wind-tunnel  curve  established  by 
Wieselsberger.  These*  incidentally*  fall  in  line  with  Shakespeare’s 
measurements  obtained  by  dropping  spheres  in  mine  shafts.  Liebster* 
reference  (9)  (0)»  performed  similar  experiments  to  those  of  Allen, 

These  data  have  relatively  high  scatter  among  themselves  and  differ 
from  Wieselsberger ' s  curve  by  as  much  as  17  percent  at  Rd  =  1.15  x 
KM.  Better  agreement  with  the  wind-tunnel  data  was  obtained  by 
Lunnon*  reference  (10)  (Q^also  by  dropping  spheres  in  water.  . 
There  is  some  discrepancy  at  the  higher  Reynolds  numbers  (Rd  >  1CH)* 
indicating  that  the  rise  in  Cd  with  Reynolds  number  is  at  a  higher 
Rd  than  in  the  wind-tunnel  tests.  It  seems  that  it  is  more  difficult 
to  obtain  a  good  drag  coefficient  in  a  water  drop  test  than  it  appears 
at  first.  In  sphere  drop  tests  at  NOL*  where  the  ca,vity  behind  the 
sphere  was  the  center  of  study  and  not  the  drag  coefficient*  it  was 
observed  that  a  sphere  never  descends  in  a  straight  line.  It  always 
falls  in  a  random  trajectory  similar  to  the  path  of  a  knuckle  ball. 
This  kind  of  trajectory  is  caused  by  the  vortices  being  shed  at 
the  separation  point  which  varies  along  the  periphery  of  the  sphere* 
thus,  altering  pressure  distribution  and  the  direction  of  the  fall. 

The  vortex  frequency  has  been  studied  by  Moller*  reference  (13) * 
and  the  Streuhal  number  correlated  with  Reynolds  number.  Excellent 
photographs  of  the  vortices  are  included  in  the  paper.  In  all  the 
early  experiments  mentioned  above*  the  length  of  fall  was  timed 
between  different  heights  without  regard  for  the  direction  of  the 
velocity  vector  of  the  center  of  gravity  at  each  point  during  the 
fall.  In  a  wind-tunnel  test*  the  drag  force  is  measured  in  the 
direction  of  the  velocity  vector;  while  in  the  drop  test  in  water* 
this  condition  was  not  satisfied.  To  obtain  an  accurate  drag 
coefficient  from  a  drop  test*  one  needs  to  determine  the  trajectory 
of  the  sphere  and  base  the  drag  calculations  on  the  actual  fall  path 
in  a  manner  similar  to  modern  data  reduction  programs  in  a  ballistics 
range . 
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Compressible  Data 


Recent  measurements  of  sphere  drag  were  performed  in  a  ballistics 
range  by  Goin  and  Lawrence,  reference  (11) *  These  tests  were  made 
in  air  and  the  compressibility  effects  investigated  by  covering  a 
Mach  number  range  between  0*2  and  0.98.  The  results  are  presented 
in  Figure  4*  In  this  figure,  Wieselsberger* s  incompressible  curve 
is  reproduced  from  Figure  3  and  the  data  at  higher  Mach  numbers 
presented  by  different  symbols*  Data  at  Mach  number  0.2  (A)  in 
a  ballistics  range  show  excellent  agreement  with  the  curve  established 
by  Wieselsberger  at  M  ^  0.1  in  a  wind  tunnel.  The  maximum  deviation 
is  at  Rd  =  3.7  x  10^  and  is  approximately  1.6  percent.  At  other 
Reynolds  numbers  where  data  were  obtained,  they  fall  directly  on 
the  Wieselsberger  curve.  The  next  set  of  data  measured  by  Goin  and 
Lawrence  is  at  M  =  0.33  (  A  ).  At  low  Reynolds  numbers  there  is  a 
very  small,  but  a  definite  effect  of  Mach  number  on  CD*  but  at 
Rd  -  104,  the  effect  is  quite  pronounced.  A  four  percent  higher 
Cd  was  measured  at  M  =  0.33  (A)  than  at  M  =  0.2  (♦).  This  leads 
one  to  conclude  that  the  compressibility  effects  are  noticeable  at 
slightly  lower  Mach  numbers  than  the  previously  accepted  value  of 
0.3*  which  was  arrived  at  by  allowing  a  change  in  density,  Ap/p, 
equal  to  0.05*  reference  (14).  The  Mach  number  range  in  these  tests 
was  extended  to  0.98.  At  M  =  0.46  (A)  and  M  =  0.60  (•)>  the  data 
cover  a  Reynolds  number  ranee  between  2  x  102  and  10^.  At  Mach 
numbers  of  0.75  (♦)>  O.89  (A)*  and  O.98  (A),  the  Reynolds  number 
coverage  is  rather  limited.  The  compressible  data  do  seem  to  follow 
the  same  slope  as  the  incompressible  curve  so  that  the  effect  of 
compressibility  (or  Mach  number)  can  be  quite  well  established  at 
least  in  a  narrow- band  of  Reynolds  numbers.  The  lowest  pressure  in 
the  ballistics  range  during  any  one  of  the  shots  was  7  torr.  It 
should  be  possible  to  measure  pressure  at  this  level  to  within  .5 
percent,  and,  with  a  three  percent  slowdown  during  the  flight,  the 
error  in  aV/dx  should  be  below  two  percent.  The  overall  accuracy 
in  the  drag  coefficient  should  be  within  +2.5  percent  at  the  most. 

There  has  been  one  more  significant  work  done  in  an  attempt 
to  determine  the  drag  of  spheres  at  subsonic  Mach  numbers.  This 
is  the  work  of  Heinrich  and  his  co-workers  at  the  University  of 
Minnesota,  references  (15)  and  (l6).  Wind-tunnel  tests  at  Mach 
numbers  between  0.59  and  0.86  have  been  performed  in  the  Reynolds 
number  range  between  103  and  2.5  x  lo3.  A  0.5-inch-diameter  sphere 
made  of  teflon  was  used  throughout  the  tests.  Data  also  were 
obtained  on  related,  oblate,  and  prolate  spheroids  at  selected  Mach 
numbers  and  Reynolds  numbers,  as  well  as  on  the  effect  of  sphere 
rotation  on  its  drag  coefficient. 

The  measured  drag  coefficients  for  spheres  have  been  included 
in  Figure  4.  Four  data  points  measured  at  Mach  number  0.39  {A) 
generally  agree,  within  the  expected  accuracy  of  measurements,  with 
the  data  obtained  in  a  ballistics  range  at  Mach  numbers  M  =  0.33 
(A)  and  0.46  (A)*  These  data  agree  also  in  the  trend  of  decreasing 
value  of  drag  coefficients  with  increasing  Reynolds  number.  As  the 
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Mach  number  increases ,  this  similarity  in  the  two  sets  of  data 
ceases.  At  Mach  number  0.66  (o),  the  wind-tunnel  data  are  below 
the  values  at  M  =  0.6  (ft)  of  the  ballistics  range  data  and, 
furthermore,  do  not  agree  with  the  trend  established  by  the  incom¬ 
pressible  measurements.  The  wind-tunnel  data  seem  to  be  independent 
of  Reynolds  number.  At  higher  Mach  numbers,  specifically,  M  =  0.74 
( 0 ),  the  data  indicate  a  reverse  trend,  i.e.,  increasing  drag 
coefficient  with  increasing  Reynolds  number.  At  still  higher  Mach 
numbers,  like  M  =  0.82  (  ^ ),  the  trend  in  the  data  seems  to  have 
reversed  again,  showing  a  slope  similar  to  that  of  the  incompressible 
curve,  although  not  as  steep;  while  at  Mach  number  0.86  (*a),  the 
trend  cannot  be  determined,  owing  to  larger  scatter  in  these  data. 

The  accuracy  of  the  wind-tunnel  measurements  mainly  depends 
on  the  accuracy  with  which  pressure  can  be  measured.  The  errors 
introduced  by  the  force  balance  used  to  measure  the  drag  force  are 
small  compared  to  the  accuracy  of  the  pressure  measurement  at  very 
low  levels.  To  obtain  the  Reynolds  number  variation  covered  by  the 
data,  the  pressure  was  the  only  parameter  varied,  since  the  same 
size  model  was  used  throughout,  and  the  Mach  number  was  kept  constant 
for  a  given  set  of  data.  The  highest  possible  error  is  at  the  low 
Reynolds  number,  since  the  pressure  there  is  very  low.  It  has  been 
estimated  by  authors  of  reference  (17)  that  the  maximum  error  at  the 
low  Reynolds  numbers  is  five  percent  and  at  the  high  Reynolds  numbers, 
one  percent.  The  scatter  in  these  data  is  as  high  as  3  to  4.5 
percent.  These  wind-tunnel  and  ballistics  range  data  of  refer¬ 
ence  (11)  are  repeated  in  Figure  5  to  better  compare  the  two  sets. 

The  wind-tunnel  data  are  shown  in  open  symbols  and  the  ballistics 
range  data,  in  solid.  Error  bars  have  been  added  to  each  data  point. 

A  linear  variation  in  the  magnitude  of  the  error  was  assumed  for 
the  wind-tunnel  data  from  five  percent  at  Rd  =  103  to  one  percent 
at  2.4  x  103.  A  constant  two  percent  error  was  applied  to  the 
ballistics  range  data.  Although  the  error  bars  in  the  two  sets  of 
data  overlap  at  certain  Reynolds  numbers,  it  is  clear  from  Figures 
4  and  5  that  the  ballistics  range  data  show  a  much  stronger  effect  of 
Reynolds  number  on  Cp  than  the  data  of  Heinrich  obtained  in  a  wind 
tunnel.  This  discrepancy  has  been  speculated  by  Heinrich,  et  al, 
to  be  due  to  tunnel  wall  interference,  model  support  interference, 
or  apparent  mass  effects,  which  are  a  result  of  a  change  in  kinetic 
energy  of  the  model  as  it  flies  and  decelerates  in  a  ballistics 
range.  These  are  certainly  valid  arguments  to  explain  the  differences 
observed  in  the  measurements,  but  it  is  doubtful  whether  they  have 
strong  enough  influences  to  cause  such  large  differences  as  are 
indicated  by  the  data.  In  a  recent  work  by  Zarin,  reference  (l8), 
the  effect  of  free-stream  turbulence  has  been  studied  on  the  drag 
coefficient  of  spheres  in  the  Reynolds  number  region  identical  to 
the  one  covered  by  the  data  of  Heinrich,  et  al,  reference  (l6), 
and  Goin  and  Lawrence,  refBrence  (11).  In  Figure  6,  one  of  the 
figures  from  reference  (l8j  is  reproduced.  It  shows  this  effect 
on  various  size  models  (the  turbulence  providing  screen  was  kept 
the  same  for  all  of  the  models,  16  x  16  x  .001  inch,  therefore, 
the  ratio  of  turbulent  eddy  size  to  model  size  was  altered  by  using 
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models  of  different  diameters) ,  The  effect  of  free-stream  turbulence 
on  sphere  drag  seems  to  be  quite  pronounced,  showing  the  reverse 
trend  in  Cj)  with  Reynolds  number  to  the  incompressible  curve  that 
was  observed  in  the  Heinrich,  et  al,  reference  (16),  wind- tunnel 
data.  Although  no  turbulence  level  measurements  are  available  for 
the  tests  in  reference  (l6),  it  is  possible  that  the  disagreement 
between  that  data  and  the  equivalent  data  obtained  in  the  ballistics 
range,  reference  (11),  is  largely  due  to  turbulence  in  the  free 
stream.  In  a  ballistics  range,  the  atmosphere  should  be  considered 
quiet  since  the  currents  and  pulsations  from  the  action  cf  the 
vacuum  pumps  used  to  evacuate  the  range  tube  to  the  desired  pressure 
level  are  small.  Furthermore,  the  countdown  procedures  usually 
require  several  minutes  between  the  time  the  pumps  are  shut  off  and 
the  time  of  fire  during  which  any  of  the  disturbances  should  damp 
out. 

Summary  of  Subsonic  Data 

In  summary,  the  subsonic  drag  data  on  spheres  have  been 
analyzed  and  the  experimental  results  compared  with  one  another. 

The  incompressible  curve  seems  to  be  rather  well  defined  by  three 
independent  methods.  At  low  Reynolds  numbers,  Allen* s  measurements 
in  water,  reference  (3),  on  air  bubbles,  and  spheres  made  of  amber 
agree  well  with  Coin  and  Lawrence *s,  reference  (11),  ballistics 
range  results.  Liebster's,  reference  (9),  drop  tests  in  water  are 
hampered  by  excessive  scatter  and  really  do  not  help  much  in 
defining  the  curve.  At  Reynolds  numbers  higher  than  800, 
Wieselsberger* s  wind-tunnel  data,  reference  (1),  agree  well  with 
Goin  and  Lawrence's  measurements  as  well  as  Lunnon’s,  reference  (10), 
water  drop  test  results.  The  latter  begin  to  deviate  from  the  wind- 
tunnel  curve  above  R^  =  2  x  10^.  Except  at  these  higher  Reynolds 
numbers,  where  data  from  drop  tests  in  water  and  air  and  the  wind- 
tunnel  experiments  differ  quite  markedly  and  ballistics  range  data 
are  not  available,  the  curve  established  by  Wieselsberger  from  his 
own  and  Allen's  measurements  in  reference  (1)  can  be  considered  to 
be  adequate  and  accurate  to  within  1.5  percent  for  defining  the 
incompressible  curve  for  spheres.  For  Reynolds  numbers  ibove 
2  x  10^,  additional  testing,  particularly  in  a  ballistics  range,  is 
necessary  to  define  the  drag  coefficient. 

The  compressibility  effects  have  been  investigated  by  two 
researchers,  each  giving  quite  different  results.  Goin  and 
Lawrence’s,  reference  (11),  measurements  follow  the  general  trend 
in  variation  with  Reynolds  number  established  by  the  incompressible 
data,  while  the  Heinrich,  et  al,  reference  (l6),  measurements  do 
fiLt.  For  this,  y*  usot  Vifl  alfl-  r.Vai  rtr-  Vail  fpf-fcp  rst-gr  'r&VH* 

in  this  case,  smaller  error  bars  than  the  ones  from  a  wind  tunnel, 
the  Goin  and  Lawrence  data  should  be  considered  to  define  the  drag 
coijffi  ciarit  at  Mach  NT  la  itaoh  nt&thvr  ahoyl-d  he 

considered  as  the  upper  limit  of  incompressibility. 
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SUPERSONIC  REGIME 


Experimental  Data 


As  it  was  in  the  case  of  subsonic  flow,  the  Reynolds  number 
regime  that  has  attracted  the  most  attention  in  previous  investi¬ 
gations  is  outside  the  Reynolds  number  region  of  interest  in  the 
falling  sphere  technique  of  atmospheric  sensing.  The  only  areas 
where  theoretical  calculations  have  been  made  of  sphere  drag  are 
in  the  very  large  Knudsen  number  region,  where  the  gas  is  highly 
rarefied  and  the  free  molecular  flow  assumptions  are  satisfied. 
Experimental  measurements  in  this  area  are  very  difficult  to  perform 
because  of  the  very  small  forces  on  the  body  and  the  necessity 
for  measuring  extremely  low  values  of  pressure.  Most  of  the  experi¬ 
mental  investigations  have  been  performed  in  the  Reynolds  number 
region  of  lo5  -  10°.  Here,  on  the  other  hand,  no  reliable  theoreti¬ 
cal  methods  have  been  developed  for  calculating  the  drag,  indicating 
the  complexity  of  the  continuum  theory.  The  range  of  interest  in  the 
falling  sphere  experiments  is  50  <  R^  <  5  x  1()3  and  Mach  numbers 
up  to  5.  This  regime  is  characterized  by  thick  boundary  layers 
and,  therefore,  large  viscous  contributions  to  the  total  drag.  The 
region  is  shown  on  a  Reynolds  -  Mach  number  plot  in  Figure  1  from 
which  one  sees  that  all  of  this  area  is  in  the  slightly  rarefied 
gas  flow  regime  of  slip  and  transitional  flows.  Slip  flow  regime 
is  defined  as  that  for  which  the  velocity  at  the  wall  is  no  longer 
zero,  as  in  the  case  of  continuum  flow,  but  retains  the  other 
properties  of  continuum  flow.  Transitional  regime  is  an  area 
between  slip  flow  and  free  molecular  flow  regimes  where  properties 
generally  are  not  well  defined  and  remain  unknown. 

Figure  7  shows  the  Reynolds  number  coverage  of  past  experiments 
available  in  open  literature  on  sphere  drag  between  Mach  numbers 
1  and  5  and  also  shows  the  range  cf  applicability  in  the  falling 
sphere  atmosphere  sensing  technique.  As  was  pointed  out  above, 
numerous  investigations  were  made  in  the  Reynolds  number  region  of 
105  -  10°.  This  is  where,  at  low  subsonic  Mach  numbers,  the 
transition  in  the  boundary  layer  of  the  sphere  took  place,  which 
was  accompanied  by  an  abrupt;  reduction  in  total  drag.  This  phenome¬ 
non  was  found  bo  diminish  at  higher  subsonic  Mach  numbers  and 
completely  disappear  at  supersonic  Mach  numbers.  This  is  illustrated 
in  Figure  8,  which  has  been  reproduced  from  reference  (8).  This 
region,  however,  is  outside  the  Reynolds  number  range  of  interest 
in  the  falling  sphere  technique  and  will  not  be  discussed  further. 
Relatively  few  investigations  have  been  made  in  the  Reynolds  number 
region  between  50  and  5  x  10^.  These  will  be  studied  and  discussed 
in  detail. 

The  most  comprehensive  investigation  covering  the  Reynolds 
number  region  of  interest  was  performed  by  Aroesty,  reference  (19). 
The  objective  of  these  wind-tunnel  measurements  was  to  gather  drag 
data  on  spheres  in  the  supersonic  rarefied  gas  regime  and  to  study 
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the  effect  of  heat  transfer  at  the  surface  of  the  sphere  on  drag, 
as  the  sphere  wall  temperature  was  reduced  from  300  degrees  K  to 
80  degrees  K  during  a  test  in  the  wind  tunnel.  A  large  number  of 
drag  measurements  were  made.  The  data  near  Mach  number  2  have  been 
replotted  in  Figure  9.  In  reference  (19)*  all  data  between  Mach 
numbers  1,619  and  2.183  were  plotted  together.  To  see  if  there  is 
a  Mach  number  effect  over  this  region,  these  3ame  data  were  sub¬ 
divided  into  smaller  Mach  number  groups  in  Figure  9»  Distinction 
was  also  made  as  to  the  number  of  wires  holding  the  sphere.  The 
main  purpose  for  such  a  plot  is  to  see  if  there  are  trends  with 
Mach  number  and  if  there  are  consistent  support  effects. 

Because  of  the  relatively  large  scatter  in  the  data,  they  were 
statistically  fitted  to  an  equation  of  the  form 


Cj)  =  A  +  +  C/R(j 


This  equation  is  represented  by  the  solid  line  in  Figure  9»  which 
is  reproduced  from  a  figure  in  reference  (19) • 

It  appears  that  the  fitted  curve  has  somewhat  more  curvature 
than  the  data  band  and  does  not  represent  the  data  in  the  Reynolds 
number  region  between  50  and  300.  This  really  is  subject  to  the 
equation  chosen  and  the  weighting  of  the  data  in  the  statistical 
method  applied  to  determine  the  constants  (A,  B,  and  C)  in  the 
equation.  The  maximum  width  of  the  scatter  band  on  percentage 
basis  is  +7  percent  from  a  faired  curve  through  the  data.  Since  the 
fitted  line  does  not  follow  the  scatter  band  directly  in  the  middle, 
the  error  bars  at  Reynolds  number  150,  for  instance,  give  a  value 

for  Cp  =  1.31  + "*"•[*  5$  and  at  Reynolds  number  700,  Cp  =  1.1 

There  are  several  reasons  given  in  reference  (19)  for  the  relatively 
large  scatter.  First  of  all,  there  is  an  uncertainty  as  to  the 
magnitude  and  effect  of  the  non-uniformities  in  the  flow  produced 
by  the  nozzle.  This  may  be  particularly  important  at  M  -  2,  where 
these  disturbances  may  affect  the  structure  of  the  flow  surrounding 
the  sphere  and  may  cause  significant  variations  in  the  base  drag. 

It  is  pointed  out  by  the  author  in  reference  (19)  that  the  nozzle 
used  for  M  -  2  test  was  the  least  uniform  of  all  the  other  nozzles 
used  for  the  higher  Mach  number  tests.  It  has  also  been  found  that 
in  low  Mach  number  tests,  when  thick  boundary  layers  are  present, 
the  results  are  sensitive  to  other  obstructions  in  the  flow  such  as 
shields  and  pressure  probes.  I-t  was  also  found  that  means  of  holding 
the  model  in  the  test  section  had  an  effect  on  the  drag  measured. 

In  all  these  tests,  the  spheres  were  wire  supported.  As  can  be 
concluded  from  Figure  9,  the  three-wire  mounted  model  consistently 
possessed  a  slightly  higher  drag  coefficient  than  the  four-wire 
mounted  model.  This  was  only  observed  for  the  smaller  size  models 
used,  while  the  one-inch-diameter  model  did  not  show  this  phenomenon. 
Another  effect  observed  was,  as  pointed  out  in  reference  (19) ,  that 
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there  was  a  non- systematic  variation  in  drag  with  the  size  of  the 
model  at  a  fixed  value  of  Reynolds  number.  These  effects  look 
very  much  like  those  studied  by  Zarin,  reference  (18),  in  the 
subsonic  regime,  where  the  drag  was  found  to  be  a  function  of 
free- stream  turbulence  and  the  relative  size  of  the  turbulent 
eddies  to  model  diameter.  Of  course,  there  was  the  ever-present 
non- uniformity  in  the  nozzle,  as  well  as  the  model  and  its  support 
assembly  effects.  All  these  contribute  to  the  widening  of  the 
scatter  band.  Except  for  a  definite  effect  of  the  number  of  support 
wires  used  on  models  of  certain  sizes,  although  overlapping  does 
occur,  no  systematic  variation  was  detected  with  either  Mach  number 
or  model  diameter.  Therefore,  the  width  of  the  scatter  uand  of 
+7  percent  should  be  considered  as  the  maximum  uncertainty  in  the 
cfata. 


Results  from  other  investigations  are  presented  in  Figure  10. 
Aroesty' s  data,  reference  (19) ,  are  represented  by  the  solid  line 
which  is  the  same  as  the  one  in  the  previous  figure  and  the  scatter 
in  these  data  represented  by  the  two  dashed  lines.  Ten  years 
before  Aroesty' s  work,  experiments  were  performed  by  Kane,  refer¬ 
ences  (20)  and  (21),  in  the  identical  Berkeley  Low  Density  Wind 
Tunnel  as  was  used  by  Aroesty.  The  reason  for  repeating  these 
experiments  later  by  Aroesty,  as  stated  in  reference  (19),  is  that 
there  was  some  uncertainty  in  interpreting  the  various  pressure 
probe  readings  in  the  low  Reynolds  number  regime.  Although  the 
force  measurements  in  the  earlier  experiments  are  reported  to  be 
quite  good,  the  values  of  Cp  and  Reynolds  number  may,  therefore, 
be  questionable.  Kane's  data  are  plotted  point-by-point  in  Figure 
10.  They  are  again  subdivided  into  smaller  Mach  number  groups  to 
check  for  Mach  number  effects.  None  were  found.  The  sphere  sizes 
varied  between  0.10  inch  to  1.00  inch  in  diameter  and  were  supported 
on  circular  section  rods  normal  to  the  flow  direction.  No  model 
size  effects  were  observed  in  these  experiments  as  well.  Kane's 
results  agree  quite  well  with  those  of  Aroesty  at  lower  Reynolds 
numbers,  where  they  were  considered  of  questionable  accuracy,  but 
show  a  systematically  higher  value  of  drag  coefficient  in  the 
higher  Reynolds  number  range  where  good  agreement  between  the  two 
sets  of  measurements  was  expected. 

Another  investigation  performed  in  the  same  Berkeley  Low 
Density  Wind  Tunnel  is  that  by  Sherman,  reference  (22).  The 
principal  reason  for  doing  these  tests  was  to  investigate  whether 
the  cross-stream  support  rod  affects  the  sphere  boundary  layer 
and,  therefore,  the  drag.  In  these  tests,  the  sphere  was  supported 
on  a  similar  rod  as  before  (references  (20)  and  (21))  but  held 
the  sphere  from  the  back,  so  that  the  rod  was  immersed  in  the 
separated  flow  region  and  in  the  wake  behind  the  sphere.  Data 
obtained  with  this  support  show  a  somewhat  lower  value  of  drag 
coefficient  at  low  Reynolds  numbers  than  the  wire-held  sphere  drag 
data  of  Aroesty  and  the  cross-stream  rod  support  of  Kane.  At 
higher  Reynolds  numbers,  Sherman' s  data  are  above  those  of  Aroesty 
and  agree  quite  well  with  those  of  Kane. 
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It  can  be  concluded  that  all  three  investigations  in  the 
identical  facility  at  Berkeley  really  show  the  same  difficulties 
in  the  measurements.  They  show  that  at  M  2  2  there  is  no  definite 
discernable  trend  with  small  changes  in  Mach  number  (1.62  <  M  <  2.8)* 
that  there  is  an  effect  of  model  size*  but  an  inconsistent  one*  and 
that  the  method  of  supporting  the  sphere  does  affect  the  measure¬ 
ments  considerably*  but  that  there  is  no  way  to  eliminate  it  from 
the  measured  data. 

There  is  one  more  reference  where  drag  measurements  on  spheres 
were  made  in  a  wind  tunnel  and  results  reported.  This  is  the  work 
of  Heinrich  and  his  co-workers  at  the  University  of  Minnesota* 
reference  (16).  These  measurements  are  shown  in  Figure  10  with  a 
square  symbol  (D*BjD)  and  are  separated  according  to  Mach  number. 
Mach  number  1.5  (O)  and  2.5  (□)  data  show  a  definite  and  consistent 
effect  of  Mach  number  on  drag*  with  M  =  1.5  having  the  higher  CD 
values.  The  Mach  number  2  (■)  data*  however*  intersect  both  the 
lower  and  higher  Reynolds  number  data.  Furthermore,  the  data  at 
low  Reynolds  numbers  are  estimated  in  reference  (16)  to  have  maximum 
error  as  high  as  27.9  percent  for  M  =  1.5  and  =  243.6  and  about 
10  percent  for  M  =  2.0  and  2.5  and  Rd  =  4-00*  thus  severely  weakening 
the  above  conclusion  of  Cp  dependency  on  Mach  number.  The  data  at 
higher  Reynolds  numbers  seem  to  agree  quite  well  with  those  of  Kane* 
references  (20)  and  (21)*  and  Sherman*  reference  (22). 

The  sphere  drag  data  described  so  far  were  obtained  in  wind- 
tunnel  facilities.  Some  tests  have  been  performed  in  a  ballistics 
range  at  the  Naval  Ordnance  Laboratory  by  May*  reference  (23) * 
and  May  and  ¥itt*  reference  (24).  In  the  former  work*  spheres  of 
sizes  varying  from  1/4  to  1/32  inch  in  diameter  were  flown.  Since 
the  smaller  size  model  could  not  be  detected  as  it  flew  in  the  range 
to  trigger  the  shadowgraph  stations*  the  spheres  were  launched  in 
clusters  that  included  some  1/4-inch-diameter  spheres  to  assure 
triggering.  In  order  that  the  spheres  would  slow  down  approximately 
at  the  same  rate,  the  ballistic  coefficient*  W/CpA*  was  matched  as 
closely  as  possible  by  using  different  density  materials  for  the 
different  size  models.  Up  to  12  spheres  were  launched  simultaneously. 
Interference  between  the  models  was  expected*  but  there  was  enough 
dispersion  between  the  spheres  that  some  of  them  were  separated  far 
enough  from  the  rest  so  they  could  be  considered  undisturbed.  Only 
the  ones  that  are  reported  to  appear  beyond  suspicion  were  replotted 
in  Figure  10.  This  was  judged  from  the  relative  position  of  the 
spheres  on  two  orthogonal  shadowgraph  plates  at  each  station  along 
the  length  of  the  flight.  Since  there  is  a  rather  large  variation 
in  Mach  number  in  these  data*  they  were  subdivided  into  two  groups 
in  Figure  10*  between  Mach  numbers  1.8l  and  2.67  (B)  and  between 
3.18  and  3.45  (8).  As  can  be  seen*  the  scatter  in  the  data  is 
too  high  to  help  to  confirm  any  of  the  wind-tunnel  measurements  in 
this  Reynolds  number  region*  where  considerable  discrepancy 
between  the  different  test  results  exists.  Neither  can  one  conclude 
from  these  ballistics  range  data  that  there  is  a  Mach  number  effect 
between  the  lowest  (1.81)  and  the  highest  (3.45)  Mach  numbers. 
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In  another  test  in  the  same  facility  by  May  and  Witt*  refer¬ 
ence  (24) ,  each  sphere  was  launched  one  at  a  time.  These  data  are 
at  higher  Reynolds  numbers  than  the  previous  set*  actually  outside 
the  range  of  interest  to  the  falling  sphere  technique  for  measuring 
the  density  of  atmosphere,  but  because  of  the  reduced  scatter  in 
these  data,  they  do  help  to  define  the  drag  coefficient  at  the 
higher  Reynolds  number  end.  As  can  be  seen  in  Figure  10,  these 
measurements  are  somewhat  lower  and  support  the  data  of  Aroesty, 
reference  (19)*  rather  than  those  of  Kane,  references  (20)  and  (21), 
and  Sherman,  reference  (22),  which  showed  higher  drag  coefficient 
values. 

Experimental  data  at  higher  supersonic  Mach  numbers  are 
presented  in  Figure  11.  At  around  Mach  number  4,  they  come  from 
two  sources,  Aroesty,  reference  (19)*  and  Wegener  and  Ashkenas, 
reference  (25).  Actual  data  points  from  the  former  reference  are 
not  shown,  but  are  represented  by  a  regression  curve  reproduced 
from  a  figure  in  reference  (19)  and  shown  by  a  solid  line,  while 
the  spread  in  measurements  is  indicated  by  the  dashed  lines.  The 
data  from  the  latter  reference  are  shown  by  an  open  circle  (O) 
which  were  also  obtained  in  a  wind  tunnel.  The  agreement  between 
the  two  sets  is  very  good.  The  Wegener  and  Ashkenas  data  are 
somewhat  below  the  fitted  curve  of  Aroesty' s  data  between  Reynolds 
numbers  1.5  x  102  and  5  x  102,  but  do  fall  almost  entirely  within 
the  lower  half  of  their  scatter  band.  Also  included  in  Figure  11 
are  the  results  of  measurements  of  around  Mach  number  6  by  Aroesty. 
The  data  are  represented  by  the  fitted  regression  curve  and  the 
scatter  in  the  data  by  dotted  lines.  The  quality  of  the  data  seems 
to  be  Improved  at  the  higher  Mach  numbers,  the  maximum  width  of 
the  scatter  band  at  M  =  6  being  about  +2  percent  from  a  faired  mean. 

Summary  of  Supersonic  Data 


Most  of  the  available  data  in  the  Mach  and  Reynolds  numbers 
range  of  interest  come  from  tests  in  wind-tunnel  facilities.  At 
lower  Reynolds  numbers,  they  come  exclusively  from  one  such  facility 
at  Berkeley.  Except  for  some  disagreement  at  the  lower  Mach  numbers 
between  Reynolds  numbers  4  x  102  and  103  with  other  wind-tunnel 
measurements,  the  data  obtained  by  Aroesty,  reference  (19)*  seem 
to  satisfactorily  represent  the  drag  coefficient  of  sphere  at 
supersonic  Mach  numbers.  In  the  region  of  discrepancy,  the  avail¬ 
able  ballistics  range  data  seem  to  support  those  of  Aroesty,  so 
that  one  may  conclude  that  Aroesty' s  data  should  be  used  throughout. 
To  increase  confidence,  it  would  be  highly  desirable  to  have  data 
generated  in  a  facility  other  than  a  wind  tunnel,  specifically,  a 
ballistics  range.  A  technique  to  obtain  drag  data  at  low  Reynolds 
numbers  in  a  ballistics  range  has  been  developed  and  demonstrated 
by  Bailey  and  Koch,  reference  (26).  This  involves  manufacturing 
models  out  of  very  low-density  materials,  in  this  cas^  foamed 
plastics  that  are  strong  enough  to  withstand  the  accelerations  in 
the  gun  during  the  launch,  causing  changes  in  the  shape  of  the  model 
or  other  of  its  physical  constants. 
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DRAG  TABLE 

Analysis  of  Existing  Drag  Table 


It  Is  Interesting  to  take  a  drag  table  that  Is  presently  used 
in  reducing  the  falling  sphere  data  to  atmospheric  densities  and 
compare  it  with  the  drag  data  that  has  been  collected  and  discussed 
in  the  previous  sections  of  this  report.  A  popular  program  in  use 
is  that  due  to  Luers  and  Engler,  reference  (27).  Figure  12  is  a 
reproduction  of  a  figure  from  reference  (27)  and  shows  a  drag  table 
generated  utilizing  Goin  and  Lawrence,  reference  (11),  in  the  sub¬ 
sonic  regime  and  Heinrich,  et  al,  reference  (l6),  in  the  super¬ 
sonic  regime.  In  Figures  13  through  15,  three  curves  representing 
high,  medium,  and  low  Reynolds  numbers  are  reproduced  from  Figure  12 
and  compared  with  the  measured  data  points  from  the  references 
discussed  earlier.  The  symbols  used  in  these  figures  are  the  same 
as  the  ones  used  earlier.  The  flagged  symbols  represent  interpolated 
values  in  cases  where  actual  data  points  in  the  specified  narrow 
Reynolds  number  range  do  not  exist,  but  can  be  interpolated  with  a 
fair  degree  of  confidence.  It  can  be  seen  that  the  drag  table 
curves  indeed  follow  the  data  of  references  (16)  and  (11)  so  that 
the  accuracy  of  the  drag  table  really  depends  on  the  accuracy  and 
validity  of  the  data  in  these  references.  As  was  pointed  out  in 
conclusions  in  the  section  on  subsonic  data,  the  Goin  and  Lawrence 
measurements  in  a  ballistics  range  are  probably  the  best  in  the 
subsonic  regime.  The  drag  table  curves  follow  these  data  very  well 
at  all  three  Reynolds  numbers  chosen  here  for  comparison  purposes. 
Supersonically,  as.  indicated  by  Luers  and  Engler,  the  Heinrich, 
et  al,  data  were  favored.  As  was  pointed  out  in  earlier  sections, 
this  probably  was  not  the  best  of  choices.  These  data  do  not 
agree  with  other  measurements,  especially  at  low  Reynolds  numbers. 
The  peaks  in  the  drag  coefficient  at  M  =  1,5  and  at  low  Reynolds 
numbers  in  Figure  12  are  solely  due  to  one  set  of  measurements, 
those  of  Heinrich,  et  al,  reference  (l6).  If  the  data  from  the 
Berkeley  Low  Density  Wind  Tunnel  are  considered,  references  (19)> 
(20),  (21),  and  (22),  then  one  would  not  conclude  that  such  a 
large  increase  in  drag  exists  at  Mach  number  approximately  1.5. 
Additional  measurements  between  Mach  numbers  1  and  2  are  highly 
desired  to  define  the  drag  coefficient  in  this  region. 

Recommended  Drag  Table 


Based  on  the  available  drag  data  in  the  open  literature,  a 
drag  table  is  recommended  as  shown  in  Figure  16.  The  drag  coeffi¬ 
cient  is  plotted  as  a  function  of  Mach  number  for  constant  values 
of  Reynolds  number,  as  was  done  in  Figure  12.  The  portion  of  the 
curves  shown  by  solid  lines  is  well  substantiated  by  experimental 
data,  while  over  the  region  shown  by  a  dashed  line,  either  measure¬ 
ments  have  not  been  made  or  the  presently  available  data  are 
contradictory.  At  subsonic  Mach  numbers,  the  data  due  to  Goin 
and  Lawrence,  reference  (11),  were  utilized  and, therefore, this 
section  of  the  plot  is  identical  to  Figure  12.  Supersonically, 
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the  Aroesty,  reference  (19),  Wegener  and  Ashkenas,  reference  (25), 
May,  reference  (23),  and  May  and  Witt,  reference  (23),  measurements 
were  considered  foremost,  and,  therefore,  the  curves  are  considerably 
different  from  those  in  Figure  12.  The  notable  difference  in  the 
presently  recommended  drag  table  is  the  absence  of  the  peak  in  drag 
coefficient  at  Mach  number  approximately  1.5.  The  measurements 
utilized  for  this  table  showed  very  little  or  no  effect  of  Mach 
number  in  the  supersonic  regime.  To  increase  its  usefulness,  the 
recommended  drag  curves  of  Figure  1 6  are  tabulated  in  Table  I. 

Recommended  Future  Tests 

In  Figure  1,  the  ranges  of  interest  in  Mach  and  Reynolds  numbers 
in  the  atmospheric  sensing  experiments  are  shown.  This  figure  is 
reproduced  in  Figure  17,  in  which  the  area  shown  by  the  line  cross 
hatching  represents  Mach  and  Reynolds  numbers  for  which  experimental 
data  are  available  and  the  drag  coefficient  of  spheres  can  be 
predicted  with  confidence.  Outside  this,  additional  tests  would  be 
desirable.  Some  incompressible  data  at  higher  Reynolds  numbers  do 
exist,  but  there  is  slight  disagreement  among  them.  It  would  be 
relatively  easy  to  extend  tests  like  those  of  Goin  and  Lawrence, 
reference  (11),  to  measure  the  drag  coefficients  there.  At  higher 
subsonic  Mach  numbers,  there  are  certain  Reynolds  number  regions 
where  interpolation  would  not  likely  give  accurate  results. 

Ballistics  range  tests  would  be  simple  and  relatively  inexpensive 
to  perform  in  this  area  and,  therefore,  are  recommended.  No  data 
exist  at  transonic  and  low  supersonic  Mach  numbers.  Tests  have 
been  made  between  Mach  numbers  1.5  and  3,  mostly  in  wind  tunnels. 
Considerable  contradiction  exists  between  different  test  results  in 
this  region  while  each  set  of  data  has  rather  large  scatter  or 
probable  error.  Aroesty' s,  reference  (19),  measurements  are  the 
most  comprehensive  and  his  data  are  recommended.  At  Mach  numbers 
around  2,  these  data  possess  approximately  +7  percent  scatter,  and, 
therefore,  it  would  be  advisable  to  confirm  these  results  with 
additional  tests  in  a  ballistics  range.  A  technique  to  obtain  drag 
data  at  low  Reynolds  numbers  has  been  developed  by  Bailey  and  Koch, 
reference  (26).  A  ballistics  range  test  has  the  advantage  of 
having  no  local  disturbances  around  the  model,  like  wire  or  string 
supports,  or  disturbances  in  the  free  stream  that  hamper  a  wind- 
tunnel  test,  especially  at  low  Reynolds  numbers,  where  they  have 
a  marked  effect  on  the  drag. 

APPLICABILITY  OF  GROUND  TEST  DATA  TO  FULL-SCALE  FLIGHT 

The  models  used  in  all  the  tests  discussed  so  far  have  been 
perfect  spheres,  or  as  nearly  perfect  as  has  been  possible  to 
manufacture  them.  The  out  of  roundness  was  checked  and  controlled 
so  that  it  would  be  within  certain  small  tolerances.  The  surface 
was  usually  polished  to  have  a  minimum  surface  roughness.  In  all 
cases,  the  spheres  were  made  of  solid  material  and  can  be  considered 
undeformable  under  the  loads  experienced  in  the  tests.  In  all  but 
one  investigation,  which  will  be  described  later  in  this  report. 
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the  surface  temperature  of  the  sphere  was  not  altered  from  that 
which  normally  exists  during  these  tests.  In  case  of  the  wind 
tunnel,  the  sphere  reaches  an  equilibrium  with  the  surrounding 
stream  and  is  considered  to  be  at  the  recovery  temperature  of 
the  free  stream.  In  a  ballistics  range,  the  model  normally  is  at 
room  temperature  when  it  is  loaded  in  the  gun  and  because  of  the 
short  flight  duration,  the  total  heat  input  during  the  flight  is 
not  high  enough  to  raise  its  temperature. 

In  an  atmospheric  sensing  experiment,  on  the: other  hand,  the 
vehicle  considered  herein  is  an  inflatable'  balloon  for  which  the 
above-mentioned  properties  may  not  apply.  These  balloons  are  made 
of  0.5-mil- thick  mylar  cut  into  panels  which  are  butt  jointed  and 
taped.  Because  of  this,  there  are  surface  irregularities  and 
roughness  elements  at  the  taped  joints.  They  may  be  deformed  due 
to  the  surrounding  pressure  distribution,  especially  at  lower 
altitudes.  Finally,  the  heat  transfer  at  the  wall  of  the  balloon 
may  not  be  the  same  as  in  a  ground  test,  due  to  solar  heating  of 
the  balloon  surface.  1  ! 

A  question  that  remains  is  how  applicable  are  the  data 
gathered  in  various  testing  facilities  on  model  spheres  to  the 
atmospheric  sensing  balloon  flights.  There  has  been  son^e  attempt 
made  to  investigate  the  effect  of  parameters  such  as  rotation, 
heat  transfer,  and  surface  roughness  on  the  drag  coefficient.  These 
will  be  discussed  in  the  following  sections  of  this  report. 

Effect  of  Rotation 

The  rocket  carrying  the  tightly  packed  mylar  balloon  6n  Its 
way  up  is  spin  stabilized.  Therefore,  it  can  be  expected  that 
after  ejection  the  inflated  sphere  will  be  rotating  as  well.  The 
spin  rate  of  the  dart  prior  to  ejection  is  between  25  and  30, rpp. 
After  the  sphere  Is  deployed  and  fully  inflated,  the  rotational 
velocity  will  decrease  considerably  because  of j the  increase  in 
moment  of  inertia.  To  estimate  the  rotational  velocity  of  the 
balloon,  let  us  apply  the  law  of  conservation  of  angular  momentum 
to  the  system  consisting  of  the  balloon  material.  It  is  assumed  ' 
here  that  the  isopentane  gas  and  the  aluminum  capsule  containing 
it  have  a  negligible  effect  on  the  moment  of  Inertia  of  the  jsystem. 
The  mylar  material  is  assumed  to  be  homogeneously  distributed  in 
a  cylinder  Inside  the  staves  in  the  initial  position  and  in  the 
walls  of  the  sphere  in  the  final  position.  Conservation  of  angular 
momentum  requires  that  , 
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where  rQ  is  the  outside  radius  and  r^  is  the  inside  radius  of  the 
sphere.1  Solving  for  Wf/u^  and  gets 
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Since  the  thickness  of  mylar  is  very  small  compared  to  the  sphere 
diameter,  all  terms  involving  *  may  be  neglected  compared  to  the 
constants  inside  the  brackets.  Therefore,- 
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There  is  only  one  work  in  which  the  effect  of  sphere  rotation  was 
specifically  measured  in  a  wind-tunnel  test.  This  was  reported  in 
reference  (16).  Figure  18  is  a  reproduction  of  a  figure  from  the 
above  reference.  Although  the  Reynolds  numbers  in  these  tests  were 
above  those  of  interest  in  the  falling  sphere  experiments,  they  may 
indicate  qualitatively  a  trend  at  lower  Reynolds  numbers  as  well. 

The  measurements  were  made  at  low  velocities,  certainly  in  the 
incompressible  regime.  In  Figure  18,  the  drag  coefficient  is  plotted 
against  the  ratio  of  the  circumferential  and  free-stream  velocities, 
«r/V.  It  can  be  seen  that  Cp  does  vary  with  the  velocity  ratio  but 
that  for  small  spin  rates,  i.e.,  for  values  of  ur/V  less  than  0.2, 
it  appeayo  tc  he"  ccnotunt .  21  is  very  unlikely  that  such  high  spin 

rates  will  ever  be  experienced  by  the  falling  balloons.  As  was  seen,  - 
the  spin  rate  resulting  from  the  dart  rotation  is  very  low. 
Furthermore,  the  fall  velocity  at  the  high  altitudes  will  be  high, 
resulting  in  negligibly  small  velocity  ratios,  wr/V.  At  low  altitudes, 
the  roughnesses  and  irregularities  in  the  balloon  surface  could 
induce  some  rotation,  but,  here  again,  the  values  of  spin  rates  will 
be  small  compared  to  the  fall  velocity.  Unless  there  is  a  drastically 
different  effect  of  rotation  on  drag  at  low  Reynolds  numbers  and 
htgEltitt  Kflnb  '  rixhi'^e  tiJtf  eil  cvr  in  FlgurA  18,  it  car  he  ccnoIuSed 

tnat  the  amount  of  rotation  that  the  falling  sphere  may  experience 
has  negligible  effect  on  its  drag  coefficient. 

Effect  of  Surface  Temperature 

It  has  been  found  that  reducing  the  surface  temperature  of  the 
sphere  in  a  wind-tunnel  test  that  is  increasing  the  heat  transfer 
at  the  wall  at  supersonic  Mach  numbers  reduces  the  drag  coefficient. 
This  has  been  reported  by  Aroesty,  reference  (19)-  In  this  inves¬ 
tigation,  the  sphere  temperature  was  reduced  to  80  degrees  K  in  a 
reservoir  above  the  edge  of  the  wind-tunnel  jet,  and  the  sphere  was 
permitted  to  fall  into  the  jet,  where  its  trajectory  was  recorded 
by  high-speed  motion  picture  camera.  As  in  a  ballistics  range  test, 
the  duration  of  the  test  is  so  short  that  it  is  assumed  that  the 
sphere  temperature  does  not  change  during  the  test. 

The  measurements  at  about  Mach  number  2  are  presented  in 
Figure  19,  which  is  a  reproduction  of  a  figure  from  reference  (19) • 
Here,  drag  coefficient  is  plotted  against  Reynolds  number  based  on 
gas  properties  behind  a  normal  shock.  Data  points  for  the  adiabatic 
wall  case  have  been  omitted;  only  the  fitted  regression  curve  has 
been  reproduced  as  a  solid  line.  This  is  the  same  curve  as  was 
shown  in  Figure  9-  The  cold  wall  data  points  are  shown  as  open 
circles  and  a  line  faired  through  these  data.  It  is  seen  that  a 
severe  c  >oling  of  the  wall  from  300  degrees  K  to  80  degrees  K 
reduces  t^e  drag  coefficient  at  the  lower  Reynolds  number  range  by 
about  nine  percent.  At  higher  Reynolds  numbers,  the  cold  wall  data 
seem  to  fair  into  the  adiabatic  data. 

In  Figures  20  and  21,  the  cold  wall  data  are  compared  to  the 
regression  curves  fitted  to  the  adiabatic  data  at  Mach  numbers  near 
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4  and  6,  respectively*  Here  again,  the  drag  coefficient  is  plotted 
against  Reynolds  number  based  on  properties  behind  the  normal  shock. 
Although  the  scatter  in  these  data  is  quite  large,  it  can  be  con¬ 
cluded  that  there  is  a  marked  decrease  in  the  drag  coefficient  as 
the  temperature  of  the  sphere  is  decreased. 

The  spheres  in  these  tests  were  subjected  to  very  severe 
cooling  conditions.  The  temperature  of  the  sphere  was  lowered  from 
room  temperature  of*  approximately  300  degrees  K  to  about  80  degrees 
K,  This  was  accompanied  by  a  reduction  in  drag  coefficient  of  as 
much  as  ten  percent  from  the  standard  incompressible  curve  values* 
During  a  balloon  flight,  on  the  other  hand,  cooling  of  its  surface 
below  the  equilibrium  temperature  cannot  occur.  In  fact,  heating 
of  the  surface  mainly  from  solar  radiation  is  expected.  There  has 
been  one  test  conducted  on  heated  spheres  at  Mach  number  near  4 
and  in  the  Reynolds  number  range  30  <  Re2  <  50,  reference  (28), 

The  sphere  temperature  was  raised  so  that  the  wall- to- stagnation 
temperature  ratio  was  2,6  and  about  ten  percent  increase  in  drag 
coefficient  resulted.  The  data  points  are  shown  by  an  open  square 
symbol  (□)  in  Figure  20,  It  is  impossible  for  the  mylar  material 
in  the  balloon  to  sustain  such  large  temperature  increases. 
Therefore,  only  a  very  small  rise  in  the  drag  coefficient  can  be 
expected  at  supersonic  Mach  numbers  as  the  balloon  gets  heated 
from  the  sun. 

The  effect  of  surface  temperature  on  drag  coefficient  at 
subsonic  Mach  numbers  and  under  incompressible  conditions  is  not 
known.  In  these  regimes,  what  governs  the  drag  coefficient  is  the 
contribution  from  the  base  drag,  which,  in  turn,  is  dependent  on 
the  location  of  flow  separation  on  the  body.  As  the  balloon  surface 
i3  heated  by  radiation  from  the  sun  and  the  surroundings,  the 
boundary  layer  may  be  stabilized  and,  therefore,  retard  separation. 
This  effect  is  not  known,  and  a  detailed  study  of  this  problem  is 
recommended. 

Effect  of  Surface  Roughness 


Unlike  the  sphere  models  used  in  the  many  tests  described 
earlier  in  this  report  that  generally  had  polished  surfaces,  the 
inflated  spherical  balloon  in  flight  will  have  some  roughness  on 
its  surface.  The  mylar  balloon  is  folded  and  squeezed  between  a 
set  of  cylindrical  staves  that  fit  inside  the  dart  used  for 
launching.  It  may  be  left  in  this  tightly  packed  configuration  for 
some  time  before  it  is  launched,  during  which  wrinkles  set  into 
the  mylar.  After  it  is  inflated,  the  mylar  may  not  be  the  smooth 
material  that  it  was  at  the  time  of  manufacture,  but  contains  tiny 
fold  marks  all  over  it.  Any  roughness  elements  will  probably  be 
small  compared  to  bound ary- layer  thickness.  Another  source  of 
roughness  is  the  seams.  Here,  the  mylar  material  is  butt  joined 
and  secured  with  heat- sensitive  tape.  This  could  be  a  significant 
source  of  roughness,  since  the  wrinkles  in  these  seams  are  quite 
pronounced.  There  also  is  some  dimpling  in  the  panels  where  the 
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radar  reflective  material  is  attached  to  the  interior  of  the  wa21. 
This  should  be  considered  more  like  a  surface  irregularity  rather 
than  roughness.  It  could  affect  the  point  of  separation  of  the 
boundary  layer  and,  therefore,  drag* 

There  have  been  some  tests  performed  on  spheres  with  surface 
roughness.  Roughened  spheres  were  launched  in  a  ballistics  range 
at  supersonic  Mach  numbers,  although  at  Reynolds  numbers  higher 
than  a  falling  balloon  would  experience*  The  results  are  reported 
by  Charters  and  Thomas  in  reference  (29)  and  are  reproduced  v.t 
Figure  22.  The  solid  line  represents  data  on  smooth  spheres 
generated  during  the  same  investigation.  The  models  used  were 
smooth  ball  bearings.  The  three  data  points  (&)  in  Figure  22 
were  obtained  with  a  rough  sphere*  In  t.hio  case,  a  9/lo-inch- 
diameter  ball  bearing  was  annealed,  and  a  criss-cross  pattern  of 
l/l6-inch-deep  grooves  were  filled  over  the  surface  about  1/8  inch 
apart.  Only  a  slight  increase,  about  three  percent,  in  drag 
coefficient  was  measured  over  the  smooth  sphere  value. 

In  the  incompressible  regime  and  at  low  Reynolds  numbers,  a 
study  of  the  effect  of  roughness  or  drag  was  performed  by  Sivier, 
reference  (30).  It  was  found  that  rather  large  roughness  elements 
have  to  be  applied  to  the  sphere  in  order  to  have  a  measurable 
effect  on  drag.  The  largest  roughness  factor,  k  (ratio  of  rough¬ 
ness  height  to  sphere  diameter),  in  those  tests  was  0.175.  Below 
Reynolds  number  500,  the  rough  sphere  drag  coefficients  were  a 
very  small  amount  below  those  for  the  smooth  sphere;  while  for 
Reynolds  number  above  500,  rough  sphere  Cd  increased  systematically 
with  increasing  roughness  factor,  k,  Very  large  increases  in  drag 
coefficient  with  roughness  were  measured  by  Seiberg,  reference  (31), 
in  Reynolds  number  region  between  600  and  1700 »  also  at  very  low 
Mach  numbers.  Three  types  of  ballc  were  tested  in  a  shccktube, 
each  representing  different  roughnesses.  Although  no  measurements 
of  k  were  made,  the  ball  materials  were  sapphire,  gunpowder,  and 
glass  and  were  classified  as  smooth,  rough,  and  very  rough, 
respectively.  The  results  from  reference  (31)  are  reproduced  in 
Figure  23.  There  is  considerable  scatter  among  the  data  which 
can  be  attributed  to  the  random  roughness  elements  and  their  random 
distribution  over  the  surface.  The  balls  may  not  have  beer,  perfect 
spheres  to  begin  with.  These  results  indicate  that,  indeed,  the 
roughness  of  the  surface  at  Reynolds  numbers  representative  of 
those  experienced  by  the  atmospheric  sensing  balloons  at  low  Mach 
numbers  affect  the  drag  coefficient  considerably.  Although  the 
size  of  the  roughness  elements  on  the  mylar  spheres  is  relatively 
small,  compared  to  the  roughness  on  Sivier's,  r-ifer;'nce  (30),  and 
Heiberg’s,  reference  (31),  models,  it  is  conceivable  that  the  seams 
and  the  irregularities  in  the  balloon  surface  could  cause  the  drag 
coefficient  to  be  slightly  different  from  the  smooth  sphere  curve 
values.  The  amount  of  correction  for  roughness  is  net  known,  and 
a  detailed  study  of  the  effect  should  be  made,  specifically  of  the 
effect  of  roughness  on  the  flow  separation  point. 
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CONCLUSIONS 

The  presently  available  drag  data  on  spheres  in  a  Reynolds- 
Mach  number  range  of  interest  to  the  falling  sphere  technique  of 
atmospheric  sensing  have  been  collected  and  analyzed.  The  following 
conclusions  can  be  drawn  from  thia  study: 

1.  The  incompressible  curve  has  been  well  established  by 
Wieselsberger,  reference  (1),  and  Allen,  reference  (3)>  and  sub¬ 
stantiated  by  recent  measurements  by  Goin  and  Lawrence,  reference 
(11).  Except  for  Rd  >  10^,  where  slight  discrepancies  exist  between 
the  results  from  the  two  sets  of  measurements,  this  curve  can  be 
considered  to  represent  the  drag  coefficient  of  spheres  for  M  <  0.2* 

2.  At  subsonic,  compressible  Mach  number's,  the  data  obtained 
in  a  ballistics  range  by  Coin  and  Lawrence,  reference  (11),  are 
recommended.  Additional  measurements  are  desired  since  there  are 
rather  wide  regimes  of  Reynolds  number  where  measurements  have  not 
been  made- 

3.  No  experimental  data  on  sphere  drag  are  available  in  the 
transonic  Mach  number  region. 

4.  Whatever  data  exist  in  the  Mach  number  region  from  one  to 
two  are  contradictory  or  have  very  large  probable  errors.  It  is 
highly  desired  to  generate  additional  data  in  this  region. 

5.  For  Mach  numbers  between  two  and  five,  CD  can  be  predicted 
with  an  adequate  degree  of  confidence.  Measurements  by  Aroesty, 
reference  (19) 1  May  and  Witt,  reference  (24),  and  Wegener  and 
Ashkenas,  reference  (25),  were  used  to  define  the  drag  curve. 

6.  Rotation  of  the  sphere  is  expected  to  be  low  enough  not 
to  have  a  measurable  effect  on  its  drag. 

7.  Large  changes  ’  n  surface  temperature  from  that  of  equilib¬ 
rium  have  produced  measurable  changes  in  Cd  at  supersonic  Mach 
numbers.  At  low  subsonic  Mach  numbers,  the  effect  of  surface 
temperature  on  drag  iz  not  known. 

8.  In  the  Reynolds  number  range  of  interest,  large  rough¬ 
nesses  on  the  surface  of  the  sphere  have  produced  large  increases 
in  drag,  small  roughnesses,  correspondingly  smaller  changes  in 
drag.  A  detailed  study  of  roughness  at  low  subsonic  Mach  numbers 
is  recommended. 
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TABLE  I 

RECOMMENDED  DRAG  TABLE  FOR  SPHERES 


Drag  Coefficients 

for  Variou 

is  Mach 

and  Reynolds  Numbers 
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M  \ 
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4000 
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.418 

.4 
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.5 
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.525 
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CO 
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• 

.8 

1.100 

.721 

.622 

.568 

.523 
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I.255 
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1.152 
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FIG.  1  RANGES  OF  INTEREST  IN  MACH  NUMBER  AND  REYNOLDS  NUMBER  IN  THE  FALLING 
SPHERE  TECHNIQUE  OF  ATMOSPHERIC  SENSING 


COEFFICIENT  OF  SPHERES  IN  INCOMPRESSIBLE,  SUBSONIC  FLOW 


FIG.  4  CRAG  COEFFICIENT  OF  SPHERES  IN  COMPRESSIBLE,  SUBSONIC  FLOW 


COMPARISON  OF  SPHERE 
IN  A  WIND  TUNNEL  AND 


PRESENT  DATA,  M  0.21 


FIG.  6  DRAG  COEFFICIENT  OF  SPHERES. IN  TURBULENT  FI 
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FIG.  8  DRAG  COEFFICIENT  OF  SPHERES  AS  A  FUNCTION  OF  REYNOLDS 
NUMBER  AND  MACH  NUMBER  IN  THE  REGION  2  x  YJ5<R  .  <9  x  105 
(NAUMANN,  REF.  (8))  d 


SUPPORTING  METHOD 
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COEFFICIENT  OF  SPHERES  NEAR  MACH  NUMBER  2,(AROESTY,  REF. (19)) 


—  1.62— -2.18  AROESTY,  REF  (9) 
ft  2.06-2.4 

t  2.4  — 2.6  KANE,  REF.  (20)AND  (21 ) 
ft  2  6—261, 

a  2.2— 2.8  SHERMAN,  REF  (22) 
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COEFFICIENT  OF  SPHERES  IN  SUPERSONIC  FLOW 


5.471  <  M  <6.039 

WEGENER  AND  ASHKENAS,  REF.  (25) 
3.8Q  <  M  <%.  325 
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COEFFICIENT  OF  SPHERES  NEAR  MACH  NUMBERS  4  AND  6 
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FIG.  12  SPHERE  DRAG  TABLE  FROM  LUERS  AND  ENGLER,  REF.  (?7) 


DRAG  COEFFICIENT  OF  SPHERES  NEAR  REYNOLDS  NUMBER  10' 


DRAG  COEFFICIENT  OF  SPHERES  NEAR  REYNOLDS  NUMBER  6  x 


COEFFICIENT  OF  SPHERES  NEAR  REYNOLDS  NUMBER 


SUMMARY  OF  DRAG  COEFFICIENT  OF  SPHERES  BASED  ON  PRESENTLY 
AVAILABLE  DATA 
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TECHNIQUE  OF  ATMOSPHERIC  SENSING 
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DRAG  COEFFICIENT  OF  GLASS,  GUN  POWDER  AND  SAPHIRE  BALLS  AT  LOW 
MACH  NUMBERS  (SELBERG,  REF.  (31)) 


